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Abstract 

We develop stock option price approximations for a model which takes both the risk of 
default and the stochastic volatility into account. We also let the intensity of defaults be 
influenced by the volatility. We show that it might be possible to infer the risk neutral 
default intensity from the stock option prices. Our option price approximation has a 
rich implied volatility surface structure and fits the data implied volatility well. Our 
calibration exercise shows that an effective hazard rate from bonds issued by a company 
can be used to explain the implied volatility skew of the option prices issued by the same 
company. We also observe that the implied yield spread that is obtained from calibrating 
all the model parameters to the option prices matches the observed yield spread. 
Key Words: Option pricing, multiscale perturbation methods, defaultable stocks, stochas- 
tic intensity of default, implied volatility skew. 

1 Introduction 

Bonds and stocks issued by the same company carry the risk of the comp any going bankrupt . 
Although this risk has been accounted for in pricing the bonds, see e.g. ISchonbucheii (120031 ) 
for a review of this literature, t here is only a l imited number of models in e quity opti o n pric - 



ing 



Carr and Linetskvl (|2006T ). ICarr and Wul (j200ot ). 



Hull et al 



(|200J) and 



Linetskvl (|200J), 



Mertonl (|1976l ). Instead, models are developed to matc h the implied volatil ity curve, which 
is observed daily: stoc hastic volatility m o dels, see e.g. iFouque et al.l (|2000l ). and jump dif- 
fusion models, see e.g. ICont and Tankovl (J200J). These models account for the default risk 
implicitly. These approaches could be improved accounting for the risk of default. When 
a company goes bankrupt and the stock of the company plunges, the out-of-the- money put 
options are suddenly very valuable. So, when these contracts are priced, the market must 
also be accounting for these crash scenarios, besides the stochastic volatility and small jump 
effects, and this is partly why the out of the money put options are worth much more than 
what the Black and Scholes model, which does not account for bankruptcy, would predict. 
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two anonymous referees and an anonymous associate editor for their valuable suggestions. 
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In addition, one might argue that the default risk should be priced consistently across 
the markets. When the likelihood of the default of a company increases (for example when 
a credit downgrade occurs), it is not unreasonable to expect that besides the increase in 
the yield spread, the trading volume in the out of the money put options increases and their 
prices go up. In fact, the more the put option is out of the money, the higher the return it will 
br ing when there is a default, so the higher the prices should be. Thi s is the under s tand ing 
of iLinetskvl (|2006l ) and the convertible bond pricing literature(see e.g. lAvache et al.l (|2003l )). 



In this note, we will use an intensity based approach to model the default risk of the 
company that issues the stock. We will also allow the volatility of the stock price to be 
stochastic. We observe that the implied yield spread that is obtained from calibrating all 
the model parameters to the option prices matches the observed yield spread (Figure 15). 
On the other hand, our calibration exercise shows that an effective hazard rate from bonds 
issued by a company can be used to explain the implied volatility skew of the option prices 
issued by the same company. 

In our model, the volatility is driven by two processes that evolve on two different time 
scales, fast and slow. The intensity of default is a function of volatility and an idiosyncratic 
stochastic component (since the credit market and stock options market can have indepen- 
dent movements) , which is driven by two other processes that evolve on fast and slow time 
scales. Our model is similar to that of lCarr and Wul (|2006l ) and the motivation for our model 



specification can be found on page 4 of that paper. The difference of our model from the Carr- 
Wu specification is that we also incorporate the market prices of volatility and intensity risks 
into our model, and we consider the intensity and volatility to be evolving on two different 
time scales. The purpose of our specification is to be able to obtain a symptotic opt i on pric e 
approximations using multiscale perturbation methods developed by IFouque et al.1 (|2003al ). 
who used this technique in the context of stochastic vol atility models. Our m odel can be 
thought of as a hybrid of the stochastic volatility model of lFouque et al.1 (|2003al ) (whose pur- 
po se is to obtain approxima t ions f or stock option prices) and the multi-scale intensity model 
of iPapageorgiou and Sircarl ([200a) (whose purpose is to price credit derivatives). But one 



should note that following the Carr-Wu specification we allow the intensity to depend on the 
volatility, and we also incorporate the market price of default risk into o ur m odel. Let us also 
point out that our default model differs from lCarr and Linetskvl (|2006l ) and ILinetskvl (|2006l ) 
since they consider a model in which the intensity of the defaults depends solely on the stock 
price. 

Our price approximation for the option prices has seven parameters (we will refer to it 
as the 7 parameter model), one of which is the average intensity, that can be calibrated 
to the observed implied volatility. We compare the performance of this model against 3 
simpler models: [5 paramete r model] a mo d el that does not account for the market price of 
intensity risk; [the model of IFouque et al.l (|2003al )] a model that does not account for the 
default risk; [3 parameter model] a model that does not account for the stochastic volatility 
effects. When we take the average intensity to be equal to the yield spread of the corporate 
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bond with smallest maturity and calibrate the rest of the parameters to the observed option 
prices we observe that the seven parameter model not only outperforms the others (which 
is to be expected) but also that it almost fits the data implied volatility perfectly for longer 
term options (9 months and over). (In our calibration we pool the options with different 
strikes and different maturities in the same pool and we have 104 data points, we consider 
m aturities from 3 month s to 2 years). The 5 parameter model always outperforms the model 
of iFouque et al.1 (|2003al ). which points to the significance of accounting for the default risk. 



Although, the three parameter model has only two free parameters (the other free parameter 
is the average intensity of default and is already fixed to be the yield spread of the shortest 
maturity bond) it fits the implied volatility surface fairly well, considering that we have 104 
data points. See Section 3.2. 

When we calibrate all parameters of each model, including the av erage intensity , to the 
option prices then all the models we consider (except the model of IFouque et al.l (|2003al ) 
since this model does not account for the default risk) fit the the data implied volatility 
curve quite well (even the 3 parameter model). But when option implied average intensity 
is plotted against the yield curve of the shortest maturity bond over the entire period of 
the data (about a year), one observes that all models except the 7 parameter model have 
unrealistically severe option implied intensities. Only the 7 parameter model almost matches 
the yield spread of the shortest maturity bond (the maturity is about two years on the last 
day of our data). Therefore, it might be possible to derive the risk neutral probability of 
default from the option prices using the 7 parameter model. See Section 3.3. Also we should 
note that the 7 and 5 parameter models have much richer implied volatility surface structure 
(see Figures [H [2 El andSJ) and Section 2.3, which explains why they fit to the option prices 
better. 

We sho uld note that we calibr ate t he parameters of our mode l to th e daily observed option 



prices as m lFouque et al.l (|2003al ) and iPapageorgiou and Sircarl (j2006l ). ICarr and Wul (120061 ) 
on the other hand perform a time series analysis to obtain the parameters of their model. 

The rest of the paper is organized as follows: In Section 2, we introduce our model 
specification. In Section [27T| we derive an approximation of the prices of the equity options, 
and in Section 2.2 we derive an approximation to the corporate bond prices, using multi- 
scale perturbation techniques. In Section 2.3 we analyze the model implied volatility surface 
structure of the models we consider. In Section 3, we calibrate the 7 parameter model and 
the other reference models to the data. 



2 A Multiscale Intensity Model 



We will use t he do ubly-stochastic Poisson process framework of 



Bremaud 



(|198ll ). see also 

Schonbucherl (|2003l ). Let (Q,7i,¥) be a filtered probability space hosting standard Brownian 



motions W = (W t , W t , W t , W t , )> whose correlation structure is given by (2.8) 
and standard Poisson process N = (N(t))t>o independent of W. These Brownian motion 
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will be used in modeling the the volatility and the intensity (see (2. 6)). We model the time 
of default, which we will denote by r, as the first time the process the time changed Poisson 
process N t = N ^ J * Xgds^j , t > 0, jumps. 

The stock price dynamics in our framework is the solution of the stochastic differential 
equation 



dS t = S t ( r t dt + a t dwl 0) -d[N t 



\ u du 



So = x, 



(2.1) 



in which we take the interest rate [rt)t>o to be deterministic as Carr and Wu (2006) and 
Linetsky (2006) do. At the time of default the stock price jumps down to zero. Note that 
the discounted stock price is a martingale under the measure IP. Let (^)t>o be the natural 
filtration of W. The price of a European call option with maturity T > is equal to 



C(t) = E 



B(t,T)(S T - K)+l {T>T} 



T t M cr{N s ;0 < s < t} 



E 



B(t,T)e-f Xsds (S T -K) + 
(2.2) 



in which B(t,T) is the price of the bond that matures at time T, whereas the price of a 
European put option is given by 



P(t) = E 



B(t,T)((K - S T )+1 {T>T} + Kl {r < T} ) 



T t Vo-{iV s ;0 < s < t} 



KB(t, T) — E 



B{t, T)e- Xsds mm(K, S T ) 



Ft 



(2.3) 



in which 5 = (St)t>o is the solution of 



dS t = St(r + X t ) + S t a t dW t {0) , S = x 



(2.4) 



From (|2,2p and (|2.4p we see that pricing a call option on a defaultable stock is equivalent to 
pricing a call optio n in a fic t itious market that does not default but whose interest rate is 
stochastic (also see Linetsky ^20061 ) ) . Using <\2.2h and (|2.3f) . the put-call parity is satisfied, 



C(t)-P(t) = S t -B(t,T)K, 



(2.5) 



as expected. 

We will use a model specification similar to that of I Carr and Wul (120061 ) and assume the 
following dynamics for the volatility at and the intensity At: 



at = a(Y t ,Z t ), \t = (ia 2 t+f{Qt,U t 



(2.6) 
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in which 




l im ^Y t )-^A(Y t ,Z t )\dt + ^dW} 1 \ Y = y, 



5c(Z t ) - VSg(Z t )T(Y t , Z t )) dt + Vdg{Z t )dW} 2) , Z 



'o = z, 



-(m - Q t ) - ^A(Q t , Ut)) dt + ^dWP, Q = q, 
6 v e J v e 

Sc(U t ) - V6g(Z t )T(Q t , U t j) dt + Vdg(U t )dW t W , U = u 



(2.7) 



Here, A and T are the market prices of volatility risk, A and T are the market prices of 
intensity risk. We will assume the following correlation structure for the Brownian motions 
in our model: 



E[wi°\w?} 
E[W1 1] M 2) ] 
E[W t {i \w! 3) ] 



Pit, 
Put, 
0, i 



1,2, 



E[W?\wl A) \ 



0,1,2, E[W t w ,Wn=0 



put, 
« w( 4 )i 



(2.8) 



0,1,2, 



for all t > 0. This model captures the covariation of the implied volatility and the spread 
of risky bonds, and it also accommodates the fact that t he credit market an d the stock 
market can show independent movements (see Section 1.1 of ICarr and Wul (|2006l ) for further 
motivation). One sho uld note that when (3 = 0, ou r model specification for the intensity 



coincides with that of 
model specification becomes that of 
of risks are bounded functions as in 
As in 



Fou que et al 



Fou que et al 



Papageorgiou an d Sircar (2006). On the other hand, when A = our 



(2003a). We assume that the market price 



(2003a|). 



Fouque et al.l (|2003al ) and iPapageorgiou and Sircarl (|2006l ) we assume that both 



a and / are bounded, smooth, and strictly positive functions. We also assume that the 
differential equation defining Z and U have unique strong solutions. We should also note 
that we think of e > 0, 5 > as small constants. Hence, th e process Y and Q a re fast 
mean reverting and Z and U evolve on a slower time scale. See lFouque et al.1 (|2003al ) for an 
exposition of multi-scale modeling in the context of stochastic volatility models. 



2.1 European Option Price 

From an application of the Feynman-Kac Theorem it follows that the price of a call option 
satisfies the following partial differential equation 



C e ' d C e ' S (t,x,y,z,q,u) =0, t<T, 

C e ' 6 (T, x, y, z, q, u) = (x - KB(t, T))+, 



(2.9) 



5 



where the operator C e,s is given by 

C e ' 5 4 -C Q + J-Ci +C 2 + V6M 1 + 5M 2 + J -Mi 
e Ve V e 

in which 



(2.10) 



f> A / s ^ 2 <9 2 _ 9 2 52 

d 2 d ~ d 

C2 ~di + 2 a2(y ' z)%2 + (/3cj2(y ' z) + f{q ' u))x fa ~ (/3cj2(y ' z) + /(<? ' n)) '' 

a d 2 d ~ d C 2 - 11 ) 

Mi = p 2 °(y, z )9i. z ) X Q^Q~ z ~ 9{z)T{y, 9(z)T(q, u) — , 



■A4 3 = pi2vV2g(z)—— + p 3 4vV2g( 



z 



dydz dqdz 

We have denoted the price of a call option by C e,<5 to emphasize the dependence on the 
parameters e and 5. 

Here, the operator £2 is the Black-Scholes operator corresponding to interest rate level 
(3a 2 (y, z)+f(q, u) and volatility a(y, z). (l/e)Co is the infinitesimal generator of the Ornstein- 
Uhlenbeck process (Y,Q). L\ contains the mixed partial derivative due to the correlation 
between the Brownian motions driving the stock price and the fast volatility factor. M\ 
contains the mixed partial derivative due to the correlation between the Brownian motions 
driving the stock price and the slow volatility factor. 5M 2 is the infinitesimal generator of 
the two-dimensional process (Z,U). M3 contains the mixed partial derivatives due to the 
correlation between the Brownian motions driving the fast and slow factors of the volatility 
and the intensity. 

Approximating the Price Using the Multiscale Perturbation Method 

We will first consider an expansion of the price C e ' S in powers of yS, 

C e ' 5 = C§ + V5Ct + bC\ + • • • , (2.12) 
and then consider the expansion of each term in (|2.12p in powers of y/e, i.e., 

CI = C 0iA; + V^C hk + eC 2)k + ■■■ , for all k E N. (2.13) 
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It can be seen from (|2.9p and (|2.12p that Cq and C{ satisfy 



Co(T,x,y,z) = (x - K) + , and 



^£ + l=d + £ 2 ) C\ = -(mi + -j=M 3 j ("/, 



(2.14) 



(2.15) 



Cl(T,x,y,z) = 0. 
Equating the 1/e and 1 / y/e terms in (12. 14ft , we obtain 

£oCo.o = 0, CqCiq + C\Cqq = 0. 



(2.16) 



Since the first equation in (|2.16j) is a homogeneous equation in y and q, we take Co,o = 
Cofl(t, x, z,u) to be indepe ndent of y an d q. (T he y and q dependent solutions have expo- 
nential growth at infinity ( Fouque et al. (|2003b )).) Similarly, from the second equation, it 
can be seen that C^o can be taken to be independent of y and q, since £160,0 = 0. When we 
equate the order one terms and use the fact that C\Ci : q = 0, we get 



-CqCxo + C2C0 — 0, 



(2.17) 



which is known as a Poisson equation (for 62,0)- (See lFouque et al.l (|2000i ).) The solvability 
condition of this equation requires that 



(£2) Co,o = 0, in which 



(2.18) 



dx 2 



dx 



(2.19) 

where (cr 2 (y, z)\, which is only a function of z, is the average (expectation) with respect 
to the invariant distribution of Y, which is Gaussian with mean m and standard deviation 
v, and (f(q,u)), which is only a function of u, is the average with respect to the invariant 
distribution of Q, which is Gaussian with mean fh and standard deviation v. All of this is 
equivalent to saying that the average (£2) of the Black-Scholes operator £2 is with respect 
to the invariant distribution of (Y, Q), whose density is given by 



1 



2ttvv 



cxp 



y — m 



v 



+ 



m 



(2.20) 



The solution to (|2.18|) is the Black-Scholes call option price when the volatility and the 
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interest rate are taken to be (a 2 (y,z)} and j3 (c 2 (y, zy) + (f(q,u)), 

0,,o(t, x, z, u) = C BS (x; vVG/^)}, /? *)) + (f(q, «)> ; T),T — t\ 

v y (2.21) 

4 x iV(di) - KB(t,T)exp (-(/?<a 2 (y,z)> + (/(g, n)))(T - i)) jV(d 2 ), 

where 

d _ \og(x/[KB(t, T)]) + (p (a 2 (y, z)) + (f(q, u)) ± | (a 2 (j/, z)))(r - g 

vV(l/,*)>(r-t) 

From (pm) and (I2TT8D . we obtain that 

C 2 , = -£q 1 (42 - <£ 2 » O),o- (2-23) 
Equating the order ^/e terms in (|2.14f) gives 

£ C 3 ,o + AC 2 , + £ 2 Ci, = 0, (2.24) 

which is again a Poisson equation for C^fi and its solvability condition requires that 

<£ 2 ) Ci, = - (AC 2 ,o) = {CiC \C 2 - (£ 2 ))> C ,o. (2.25) 

Observe that we have Ci^T, x, z, u) = 0. We can compute Ci,o explicitly. To facilitate 
this calculation we will need to introduce two new functions. Let (f>(y, z) and (p(q, u) be the 
solutions of 

C (j)(y,z,q,u) = a 2 {y,z) - (a 2 (y, z)) and C <j>(y, z,q,u) = f(q,u) - (/(?,«)) , (2.26) 
respectively. Now, one can easily deduce that 

C \C 2 - <£ 2 »C ,o = • Q* 2 ^ + " PCofl) + * ■ " 0>,o) • (2.27) 

On the other hand, using it can be seen that 

Ci,o = (■j=vp l (a</>j\ x ^ x2 -^r) + P(^f>i(°4>t)) x2 ^-^r + (^»« W»>) l2 ^r 

- WS <A, S ) (ix^ + „Jg> - ^) - fl V2 (A*,) - C , ^ 

(2.28) 

in which and <p v are the derivatives of (f> and 0, respectively, with respect to the y variable. 
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The solution of (|2.28[) can be explicitly computed as 



Ci,o = -(r-t)RHS, (2.29) 

in which RHS is the right-hand-side of (|2.28p . using the fact that (£2) and x n d n /dx n commute 
for n G N. 

Next, we consider f|2. 15j) together with (|2.13p . Equating 1/e and l/\/e order terms gives 

£ C ,i = 0, £ Ci,i = 0. (2.30) 

The first equation in (|2.30p implies that Co 1 is independent of y and q, and the second 
equation can be obtained by observing that M^Cqa = £iCo,i = 0, which implies that C\ t \ 
is independent of y and q as well. 

Matching the first order terms gives 

C C 2A + £ 2 Q),i = -MiCofl, (2.31) 

where we have used the fact that MzCi$ = J0%C% t i = 0. Equation (|2.3ip is a Poisson equation 
and its solvability condition gives 

(£2) C ,i = - (Mi) C , , C 0j i(T,x,y,z)=0. (2.32) 

Using the facts that 



|co,o = (T - t) (a\y, z)\ Q* 2 ^ + 0x^C o>o - /3C ,o 
^C , = (T-t)(f( q ,u))Jx^-C , 



(2.33) 



in which 

(a 2 (y, z)) z = §~ z <^ 2 (y, z)) , (f(q, u)) u = |- (f(q, u)) u , (2.34) 
and the definition of the operator M 1 , we can derive the solution to (|2.32p as 

Co,i(t, x, y, z) = (T 2 - p 2 2)) <f7 2 (y, ^)> 2 (x^ (* 2 ^-^) + f 3 ^^^) 



— o — (/>« • 5 • (r) • x—— - c ,o 



2 

dx 

(2.35) 

again using the fact that (£2) and x n d n /dx n commute for n £ N. 
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Our approximation to the price of the call option then becomes 



C<,*A Coi0 _ (r _ ()(F . ^ 



d 



dx 2 



+ Vn { X 



dx 



0.0 



d 2 Q),o 
dx 2 



dC 0fi 
dx 



Co,o 



in which 



V4 



UP! (<T0 y ) , 



w(A<£, 



u <A0 y ) 



2e ( v (A0y) /? + w ^A^ + t5 ^A^ /3 + 5 ^A0 5 




2V 1 s = VSp 2 (a)g(a 2 ) z , 

2Vi = V6 P2 .(a).g. (a 2 ) z • fi - v ^^k 

2^ = • (r) • (<7 2 ) 2 • /? - (/)„ • 5 • (f ) . 



(2.36) 



(2.37) 



When we calibrate (daily) to the observed option prices (12. 36ft will be our starting point, and 
we will determine the 7 unknowns, 



X(z,u)^p{a 2 (y,z)) + (f(q,u)), 



(2.38) 



which appears in C fi , V{, V 2 e , Vi, Vf, Vi and Vi from the observed option prices. Alterna- 
tively, we will obtain A from the yield spread of defaultable bonds (either by using the yield 
spread of the shortest maturity bond directly or by calibrating the model implied bond price 
to the yield spread curve) and calibrate rest of the parameters to the observed option prices. 
In either case, the average volatility term in (|2.2ip will be calculated from the observed stock 
prices. 



Not e that when (3 = and / = we obtain the approx imation formula in iFouque et al 



(|2003al ) as expected. Moreover, as in IFouque et al.l (|2003al ). the accuracy of our approxima- 
tion can be shown to be 

\C^ S - CM| < C(e\ log(e)| + 5), (2.39) 

for a constant C > 0. The proof can be first done for smooth pay-off functions by using 
a higher order approximation for C e ' S and deriving an the expression for the residual using 
the Feynman-Kac theorem. The proof for the call option and other non-sm o oth pa y-offs can 
be performed by generalizing the regularization argument in IFouque et al.l (|2003bl ). This is 
where we use the boundedness assumption on the functions / and a. 



Remark 2.1 The approximation for the put option can be similarly derived. Note that the 
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put-call parity relationship holds between the approximate call price and the approximate put 
price. 



Note: 

In the calibration exercises below, we will compar e the pe rformance of [the seven param- 
eter model] in (|2.36[) with [the model of lFouque et al.l (|2003al )]. which can be obtained in our 
framework by setting the intensity A in (|2.6|) to be zero. In this case stochastic volatility is 
the only source of the implied volatility smile/smirk. We will also consider the performance 
of two other models: [The 5 parameter model] is obtained when we set the market risk func- 
tions to be equal to zero, i.e., A = A = r = T = 0. In this case the approximation formula 
becomes 



C*> s ± C ,o — (T — t) 
+ (T — tf 



( 2 d 2 C ,o 
dx 2 



V ^Tx < X 



+ Vfx 



2^ 2 Co,0 



dx 2 







yO x ^_ | x 2 ® C fl \ y6 x 2 

1 dx \ dx 2 J 2 dx 2 



(2.40) 



and the constants A (the average default intensity that appears in the formula for C^o), 0, 
Vf, V£ an d yf are to be calibrated to the observed option prices. 

[The three parameter model] is obtained when we set A = A = r = T = and also take 
at = o~, for some positive constant a. In this case the approximation formula becomes 



CM 4 C ,o + (— (T - t)V e + (T - tfV 



6\ j2 d 2 Co,o 
dx 2 



(2.41) 



and the constants A (the average default intensity that appears in the formula for Ceo), V e , 
V s , which are, again, to be calibrated to the observed option prices. In this last model the 
implied volatility skew is driven only by the default intensity. 

Instead of calibrating A to the option prices along with the other parameters, one might 
prefer to obtain it from the yield spread of the corporate bond prices. 



2.2 Bond Price 

Under the zero recovery assumption, the defaultable bond price is given by 



B e ' (t, Y t , Z t , Q t , U t ; T) = B(t, T)E 



cxp 



T 



((3a 2 (Y s ,Z s ) + f(Q s ,U s ))dt 



Yt, Zt, Qt, Ut 
(2.42) 

Here, we assume the interest rate to be constant because we are interested in bonds with 
short maturities. The bond price satisfies 



B e > s (t, x, y, z, q, u) =: B(t, T)B € ' 5 (t, x, y, z, q, u) 



(2.43) 
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in which B e ' satisfies 



C e ' S B e ' 5 (t,x,y,z,q,u) = 0, t < T, 
B £ > S (T,x,y,z,q,u) = 1, 



where the operator £ e ' is given by 



^ i i ~ 

C e > 6 = -C + -pCi + C 2 + V6Mi + 8M 2 + \ -Mi 
e y/e V e 

in which Co, Mi and M3 are as in (|2.1ip and 



(2.44) 



(2.45) 



Ci 4 -A( y , zjw^ - A(g, n)5V2^-, 

C 2 ±^- t -{(3a\y,z) + f( q ,u))., (2-46) 

9-9 

-A*i = -9{z)T(y,z)— - g(z)T(q,u) — , 

Let us obtain an expansion of the price B e,s in powers of y/5, 

B € ' s = Bq + VSB{ + &B| + • • • . (2.47) 
We will then consider the expansion of each term in (|2.47ft in powers of \/e, i.e., 

B% = B ,k + VZBi,k + e^2,fc + • • • , for all k € N. (2.48) 
The function Bo t o(t, z, u) can be determined from 

'£ 2 \ 5 ,o = 0, B 0fi (T, z, u) = 1, (2.49) 



as 

Bo,o(*, 2, u) = exp(-A>, u){T - t)), (2.50) 
in which A is as in ()2.38p . On the other hand, Bi$ can be determined using 

l Z 2 )Bi fl = (CiQ l ( K C2-(c 2 )))Bofi, Bi, (T,z,u) = 0, (2.51) 
the solution of which is 

Bi fi (t, z, u) = L(z, u)(T - t) exp (— A(T - t)) , (2.52) 

in which 

L(z, u) = vV2 (p {4> y K) + (^A)) + vV2 [p (4> q l) + (^X)) . (2.53) 
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The function Bq \ can be determined from 



(C 2 ) B 0> i = - (Mi)Bo fi , B , 1 (T,x,y,z)=0. (2.54) 
The solution of this equation can be written as 

B 0A {t, z, u) = -L{z, u) (r ~ t)2 exp(-A(z, u){T - t)), (2.55) 

in which 

L(z,u)=g(T){~X) z +g(r){X) q . (2.56) 

It can be shown that 

B e ' 5 (t, z, u; T) = B{t, T) exp (\ ■ (T - t)) (l + L ■ (T - t) - L ■ ^ ~ Q ^ , (2.57) 



satisfi es j.B 6 '' 5 - 5^1 < C(e + S), for some constant C > 0, as in iPapageorgiou and Sircarl 



(2006). Since very short term maturity options are not available, we will not use this formula 
to calibrate A from the yield spread data. Instead, in Section 3.1, we will take A to be equal 
to the yield spread of the shortest maturity bond. 

2.3 The Approximate Implied Volatility 

The implied volatility of the approximate pricing formula (|2.36p . which we will denote by /, 
is implicitly defined as 

C BS (x; I, r-K,T-t) = C e > 5 {t, x, z). (2.58) 

In order to understand how a typical implied volatility surface of our model looks like, we 
expand the implied volatility / in terms of e and 5 as 

I = h + sfel{ + VSI( + ■■■ (2.59) 

We then apply the Taylor expansion formula to the left-hand-side of (I2.58D and use the price 
formula (|2.36[) on the right-hand-side and obtain 

C BS (x; Jo, r; K,T-t) + (yfeli + Vdlf h,r; K, T - t) + ■ ■ ■ 

_ Co , -{T — t) (^| (**fM) + v{x ^ + y, _ Co ,„)) (M0) 

+ (T - t f (V?*l (x^) + + Vi (Jg - C, )) + ■ ■ ■ . 

Matching the order one terms, we get Iq as the solution of 

C BS (x; Jo, 0; KB(t, T),T -t) = C BS (x; y/&*j, A; KB(t, T),T - t) . (2.61) 
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Note that Iq is already a non-linear function of the log moneyness, i.e. \og{K/x), and of the 
time to maturity, i.e. T — t. In fact, 

xy/T=tN'{d x )-Q± = B(t, T)(N(d 2 ) - N(d 2 )) < 0, (2.62) 

i.e. Iq is a strictly decreasing function of K. Here, 

, A ^g{x/[KB{t,T)})±\ll-{T-t) 

The first term in the implied volatility expansion, Iq, already carries the features of a typical 
smile curve and captures the fact that in the money call contracts or out of the money put 
contracts have larger implied volatilities. 

To match a given implied volatility curve, one will need to use yS and y/e terms in the 
expansion of the implied volatility. Each of these terms will have one free parameter that can 
be used to calibrate the model to a given implied volatility curve. Matching the \fk and y/5 
terms in (j2.60p . we obtain 

rre rrr A {rre 9 f 2 dC , \ 2 ^Cofl f dC 0fi \\ f OC^y 1 a 

^ = " (T - ^V^Tx V -a* J + v > x -tor + v * I - c 0fi \ ( -^r 1 and 



Vsif - cr - t f (vf.± (^) + ^ + vi - Co,) ^ ^ 



-1 



(2.64) 



To derive these terms more explicitly we will use the derivatives of the Black-Scholes function 
Cbs on the left hand side of equation f)2 .61|) : 

~~do~ = ( *Kox -g^g- (2.65) 
^)^zC BS =(l--^=)—C BS and, (2.66) 



dxjdcr V IoVT^tJ Oct 

dC BS _ xexp(-4/2)y/T=t . . 

The corresponding derivatives of Co,o can be obtained by replacing d\ by d\ and Iq by \J (cr 2 )- 
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Now, we can obtain the correction terms in the implied volatilities as 




di 



Vfy/2*. (T-t) exp ( ^ 



vV)-(r-t). 

-B(t,T)N{d 2 ) 



rre 1 ( ~di + d\ 

V 2 -j-^ exp 1 



y/ll{ = (T-t)Vfexp 



-d~x + d\ 



+ (T — t)Vgy/2ir -(T-t) 




B(t,T)N(d 2 



<(T») 



+ (T-t)V 2 5 — exp 



r 2 



(2.68) 



Typical implied volatility su rfaces corresp o nding to the models/approximate prices in 
(|2.41|) : (|2.40p ; and the model of iFouque et al.l (|2003al ). which is the same as setting A = 
in (|2,40p ; and (|2.36p are given respectively by Figures [Q [H [3] and Although all of the 
models are able to produce the well-known properties of the implied volatility surface of the 
observed option prices (a-the skewedness (out of the money put options have more implied 
volatility then in the money ones) and b- decrease in the degree of skewedness as the time to 
maturity grows, c-smile (the minimum of the implied volatility is at the forward stock price 
SfB(t,T))) to some degree, the models in (|2.40p and (|2.36p are able to produce surfaces with 
richer structures that looks more like the real implied volatility surfaces. 



3 Calibration to Data 



3.1 Data Description 

We calibrated (daily) our model's implied volatility Ford's options over the period of 2/2/2005- 
1/9/2006. The sources of data are: 

• The stock price data is obtained from finance.yahoo.com. 

• The stock opti on data is from Opt ionMetrics under WRDS database, which is the 



same database ICarr and Wul (|2006l ) use. It is contained in their Volatility Surface 
file, which contains the interpolated volatility surface using a methodology based on 
kernel smoothing algorithm. The file contains information on standardized options, 
with expirations 91, 122, 152, 182, 273, 365, 547 and 730 calender days and different 
strikes. There are 130 different call (put) options on each day. Exchange traded options 
on individual stocks are American-style and therefore the price reflects an early exercise 
premium. Option metrics uses binomial tree to back out the option implied volatility 
th at explicitly ac c ounts for this early exercise premium. We adopt the practice employed 
by lCarr and Wul (|2006l ) for estimating the market prices of European options: We take 



15 



the average of two implied volatilities at each strike and convert them into out-of-the- 
money European option prices using the Black-Scholes formula. 



• The treasury yield curve data for each day is also obtained from OptionMetrics. It 
is derived from the LIBOR rates and settlement prices of CME Eurodollar futures. 
(Recall that we need the bond price B(t,T) for short time to maturities T — t.) 

• Ford's bond data is obtained from a Bloomberg terminal. 

We will consider the performance of the three models given in (|2,36[) , (|2,4U[) and (|2.4ip , in 
fitting the implied volatility curve of the observed option data. In the price formula (|2.36p , we 
have seven parameters to calibrate: A, Vf, Vf, i S {1, 2, 3}; on the other hand (|2.40p we have 
five parameters to calibrate: A, Vf, Vf, i G {1, 2}; and finally, in (|2.4ip we have 3 parameters 
to calibrate: A, V e , V s . We will also consider the approxim ation formula when w hen we set 



A = (this is the approximation is the same as the one in iFouque et al.l (|2003al )). and test 
the performance of the other models against it. The average volatility term (a 2 (y,z)') term 
that appear i n our ap proximation formula will be estimated using the stock price data as in 



Fouaue et al.l (|2003ah 



Following the suggestion of lCont and Tankovl ([2004]) (see page 439), we perform the cali- 



bration on a particular day by minimizing 

2 



£ (j^™ 1 ^, Ki) - I modcl (Ti, Ki- model parameters) 

i 

(O obscrvcd (Ti, Ki) - O modcl (Ti,Ki; model parameters))' 



(3.1) 



E 



veg& 2 (Ti,Ki 



in which O observed (Tj, .fQ) is the market price of the out of the money option (put or call) 
with strike price Ki, and maturity Tj, where as O model (Tj, Kf, model parameters) is the model 
implied price of the same option contract. The term vega(Tj, Ki) is computed using (|2.67p 
with Io(Ti,Ki) as the market implied volatility 

We will group our calibration scheme into two depending on whether the average hazard 
rate A, defined in (12.38H (which enters into the approximate price formulas through Cb,o in 
(I2.2ip is inferred from the yield spread of the bond prices or it is inferred as a result of 
calibration to the market option prices. These two groups of calibration schemes will be the 
topic of the following two subsections. 



3.2 When the Hazard Rate A is Inferred Prom the Bond Yield Spread 

We fix the average hazard rate A that appears in our price approximation formulas (|2.36|) . 
(|2.40p and (I2.4ip . to be the yield spread of the shortest maturity corporate bond (of the same 
company), which is 2 years in our calibration exercise. We pool the options of maturities 
91, 122, 152, 182, 273, 365, 547, 730 together and use ([53} to calibrate our model. Each matu- 
rity has 13 strike prices available, therefore our data set 104 observed option prices available. 
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After we calibrate our models and the model of lFouque et al.l (|2003al ). which is obtained by 
setting A = in (|2.40p . we plot their implied volatilities against the implied volatilities of the 
data. These are illustrated in Figures El and El Each figure is for a different maturity. 

As expected the 7 parameter model of (|2.36p fits the data implied volatility better. (Ob- 
serve that since we fixed A to be the yield spread of the shortest maturity bond there are 
only 6 parameters to calibrate). But the 7 parameter model does not only outperform the 
other models but also fits the data implied volatility almost perfectly. 7 parameter model 
outperforms the other models by a lot for very small and very large strikes. Another inter- 
esting observation is that the 5 parameter model of (|2.40p always outperform s (again since 



the A i s set to be the yield spread of the shortest maturity bond) the model of lFouque et al 



(|2003al ) which is the same as (|2.40p except that the A is set to zero instead of being set to the 
yield spread of the shortest maturity bond. This shows the significance of taking the yield 
spread of the corporate bonds into account in pricing options. Note that although the model 
of (|2.4ip had only two free parameters, it performed fairly well (considering that we had 104 
data points available). Also note that, the performance of the 7 parameter model deteriorates 
as the time to maturity becomes smaller. This is because for the shorter maturity options 
one should use shorter term yield spreads. 



3.3 When the Hazard Rate A is Inferred from Market Option Prices Di- 
rectly 

In this section, we will analyze whether it would be possible to infer the risk neural expected 
time of default if we calibrate the models of (I2.36p . (12,401) and (|2.4ip to the data implied 
volatility. That is along with the other parameters, if we also let A be decided by the prices 
of the options would it be possible to say something about the yield spread of the corporate 
bonds. 

Since each model has one extra free variable, namely A, to fit to the data implied volatility, 
not surprisingly, the performance of each model enhances, see Figures [TOl [TT1 and [T2l But 
the enhancement is not an ordinary one: all of our three models fit the data exceptionally 
well, even the three parameter model in (I2.4ip . The improvement in the most advanced 
model is the smallest. But then when we compare the data implied yield spread versus 
the yield spread of the shortest maturity bond, see Figures [131 [H] and [151 the story is a 
little bit different. The implied yield spread, A of the 3 parameter model of (|2.4ip is much 
larger than that of the yield spread of the shortest maturity corporate bond. Since, the 3 
parameter model does not account for volatility, the implied yield spread is drastically severe. 
For the 5 parameter model the difference seems to be less severe, since the model accounts 
for the stochastic volatility effects, but the difference between the yield spread curves is still 
noticeable. On the other hand, the implied yield spread for the 7 parameter model of (|2.36p 
seems to fit the yield spread of the corporate bond almost perfectly. The difference between 
the 7 parameter and the 5 parameter values is that the former one accounts for the market 
price of intensity risk. 
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Figure 1: The Implied Volatility Surface of the 3 parameter model in (|2.4ip . Average volatil- 
ity^. 2, average intensity=A = 0.02, interest rate(r)=0.04 (assuming B(t,T) = e r ( T_s )) 
V e = 0.0015, and V s = 0.001 



0.8 -, 
O.Y - 




Figure 2: The Implied Volatility Surface corresponding to the 5 parameter model in (|2.40p 
Average volatility =0.2, A = 0.02, r=0.04, V{ = -0.0015 V£ = 0.001, Vf = -0.001 V s = 
-0.001. 
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Figure 3: The Implied Volatility Surface corresponding to the model of IFouque et al. (|2003ah . 
which is the same as (|2.40p when A = 0. Average volatility =0.2, r=0.04, V{ = -0.0015 



V? = 0.001. 



-0.001 v d 



-0.001. 




time to maturity 

moneyness 



Figure 4: The Implied Volatility Surface corresponding to the 7 parameter model in ()2.36|) . 
Average volatility =0.2, A = 0.02, r=0.04, V{ = -0.0015 V 2 e = 0.001, V{ = -0.005, Vf = 
-0.001 V s = -0.001 and Kf = -0.06. 
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Figure 5: January 9th, 2006; Maturity= 273 calender days (9 months), Stock price=8.62, 
historical volatility: 29.22%, 13 strikes, interpolated LIBOR rate=4.7771%, A = 4.385% is 
taken to be the yield spread of the shortest maturity bond, which matures on which matures 
on 1/15/2008 (726 days). 
Legend: 

'o': observed implied volatility; 'x': 3 parameter model implied volatility in (12.41D : '*': 5 
parameter model im plied volatility in (12.401) ; full circle: 7 paramater model implied volatility 



in (|2.36p : diamond: iFouque et al 
A = in (12^01) . 



)'s implied volatility, which is obtained by setting 



21 



0.7 



0.65 



0.6 



■K °- 55 <b 



0.5" 



0.45 



0.4 



8 



9 



o 



0.35 



20 



30 



40 50 60 

put delta 



70 



80 



Figure 6: January 9th, 2006; Maturity= 365 days, Stock price=8.62, historical volatility: 
29.22%, interpolated LIBOR rate=4.7795%, 13 strikes, A = 4.385% is taken to be the yield 
spread of the shortest maturity bond, which matures on which matures on 1/15/2008 (726 
days). The legend is the same as that of Fgure [5j 
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Figure 7: January 9th, 2006; Maturity= 547 days (1.5 years), Stock price=8.62, historical 
volatility: 29.22%, interpolated LIBOR rate=4.7550%, 13 strikes, A = 4.385% is taken to be 
the yield spread of the shortest maturity bond, which matures on which matures on 1/15/2008 
(726 days). The legend is the same as that of Fgure [H 
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Figure 8: January 9th, 2006; Maturity= 730 days (2 years), Stock price=8.62, historical 
volatility: 29.22%, interpolated LIBOR rate=4.7357%, 13 strikes, A = 4.385% is taken to be 
the yield spread of the shortest maturity bond, which matures on which matures on 1/15/2008 
(726 days). The legend is the same as that of Fgure [5j 
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Figure 9: January 9th, 2006; Maturity= 273 calender days (9 months), Stock price=8.62, 
historical volatility: 29.22%, 13 strikes, interpolated LIBOR rate=4.7771The legend is the 
same as that of Fgure [5j 
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Figure 10: January 9th, 2006; Maturity= 365 days, Stock price=8.62, historical volatility: 
29.22%, 13 strikes, interpolated LIBOR rate=4.7771%. The legend is the same as that 
of Fgure [5j 
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Figure 11: January 9th, 2006; Maturity= 547 calender days (1.5 years), Stock price=8.62, 
historical volatility: 29.22%, 13 strikes, interpolated LIBOR rate=4.7771%The legend is 
the same as that of Fgure [5j 
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Figure 12: January 9th, 2006; Maturity= 730 calender days (2 years), Stock price=8.62, 
historical volatility: 29.22%, 13 strikes, interpolated LIBOR rate=4.7771%. 
'o': observed implied volatility; 'x': 3 parameter model implied volatility in (12.411) : '*': 5 
parameter model im plied volatility in (12.401) ; full circle: 7 paramater model implied volatility 

2003al )'s implied volatility, which is obtained by setting 



in (|2.36p : diamond: 
A = in (|2^0D . 
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Figure 13: The implied yield spread when we calibrate the 3 parameter model (I2.4ip to the 
option data with maturities of 152, 182, 273, 365, 547, 730 (each maturity has 13 strikes) 
days versus versus the yield spread of the shortest maturity corporate bond. 
Option implied A is solid blue line, yield spread is green broken line. 
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Figure 14: The implied yield spread when we calibrate the 5 parameter model (I2.40P to the 
option data with maturities of 152, 182, 273, 365, 547, 730 (each maturity has 13 strikes) 
days versus versus the yield spread of the shortest maturity corporate bond. 
Option implied A is solid blue line, yield spread is green broken line. 
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Figure 15: The implied yield spread when we calibrate the 7 parameter model (|2.36|) to the 
option data with maturities of 152, 182, 273, 365, 547, 730 (each maturity has 13 strikes) 
days versus versus the yield spread of the shortest maturity corporate bond. 
Option implied A is solid blue line, yield spread is green broken line. 
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